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The  Vibrational  Relaxation  Processes  in  a  CO  -N  -H  0  Laser 
System  .  222 


by  Yan  Haixing 

(Institute  of  Mechanics,  Academia  Sinica) 

Paper  received  op^November  6 ,  1979/  v,  ,Cf  Jr' 

y  <v" 

Vibp^fiqnal  relaxation  processes  and  data 
of  a  f(X§-K£-H&  laser  system  up  to  the  present 
are  suTvfiyeo.'  The  data  before  1988  are  based  on 
Taylor  &  Bitterman's  review.  The  data  are 
treated  with  weighted  least  square  method  to  get 
the  best  fit  in  the  form  of  log^X*— A+BX+OX*+DX* 
and  the  corresponding  coefficients  are  ^  ° 

presented.  Using  these  best  fits  of  data  and 
relaxation  model  with  three  modes  and  four 
vibrational  temperatures  which  is  based  on  the 
present  work,  numerical  calculation  results  of 
small  signal  gain  of  GDL  agree  with  experiments, 
thus  eliminating  the  error  introduced  by  the 
uncertainties  of  relaxation  data  and  model.  /«-•*> 

-f  >  rtt'O  1  r<-  i  f  1  - 

I.  Preface  ehj 


The  relaxation  process  of  the  vibrational  degree-of-f reedom 
inside  molecules  has  a  great  effect  on  the  fields  of  high  speed 
gasdynamics  [33,  atmospheric  physics  £43,  acoustics  [53, 
chemodynamics  (6~83  and  so  forth.  In  the  past  ten  years,  a  series 
of  research  concerning  the  vibrational  relaxation  process  was 
widely  conducted.  A  series  of  papers  [9  ~  131  as  well  as  review 
articles  Cl,  14  ~173  were  presented.  The  relaxation  process  has 
created  a  key  effect  on  gas  laser  application.  Without  sufficient 


knowledge  about  this  relaxation  process,  it  is  impossible  to 
invent  a  new  laser  system,  to  improve  device  performance  or  to 
conduct  analysis  of  the  laser  dynamic  process.  Laser  technology 
development  requires  in-depth  understanding  about  the  relaxation 
process;  however,  it  itself  has  become  a  powerful  measuring  tool 
for  the  relaxation  process  experiment. 

At  the  beginning  of  1969,  a  review  concerning  the  vibrational 

relaxation  process  of  a  CO  -N  laser  system  was  presented  by 

2  2 

Taylor  and  Bitterman  Cll.  This  review  is  still  widely  recommended 

today;  however,  it  lacks  data  for  many  important  relaxation 

processes  because  of  the  research  level  at  that  time.  Data 

process  and  analysis  were  imperfect  which  made  some  of.  the  rate 

data  inaccurate;  moveover,  some  important  problems  in  the 

relaxation  process  of  the  system  were  not  thoroughly  studied. 

Since  1969,  a  wealth  of  data  were  accumulated  for  the  relaxation 

process  of  the  CO  -N  -H  0  laser  system.  In  order  to 

2  2  2 


calculate  the  performance  of  a  gas  dynamics  laser  CGDL),  some  new 
data  were  collected  and  used  [18,191;  however,  no  detailed 
analysis  for  the  relaxation  process  was  performed.  The  quantity 
of  collected  data  was  also  little;  only  some  very  imperfect  rate 
values  were  presented,  and  thus  the  calculation  results  of  the  GDL 
performance  showed  large  errors  C 21 . 


Based  on  Reference  [11,  all  new  experiment  data  and  analytical 
results  from  then  until  the  present  are  further  collected  and 


processed;  the  best  curve  fits  and  analytical  formula  for  various 

'■  s' 
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C  i<\.X 
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relaxation  rates  are  presented,  and  thus  can  be  easily  to  used. 

The  best  curve  fits  presented  in  this  paper  extend  the  data 
applicable  range,  and  point  out  the  major  difference  from  the  Ref. 
til  results.  Good  results  are  also  obtained  by  applying  this 
paper's  results  to  the  non- equilibrium  calculation  of  a  burning 
type  GDL.  This  means  that  the  reliability  of  this  paper's  results 
is  again  proved  from  an  additional  aspect. 

II.  Relaxation  Mechanism 

The  energy  exchanging  reactions  considered  here  are  Csee  Fig. 

1)  : 


V-T  exchange: 


1/cm  Cl) 


OKfc) + v  ~tsM  av-hSt+«e7  l/cm 


Nt+HfO 


■  lft+HjO+2881  1/cm 


C  2) 


1/cm  (3) 


1/cm  C  4) 


1/cm  C5) 


C  6) 


Intermolecule  V-V  exchange: 


005(r*).+N*  ♦=i;CO,+N;4-18  1/cm  (7) 

V-V  exchange  inside  molecule: 

0&5(«*)+<Xfc  3EapDC<i%*f»0+OQi+aW  1/cm  (8) 

a>J(*,)+Nt  ZliOD5<i^+*t)+N,+272  1/cm  (9) 


COSO*) 4W  5=S00Kn+*p) l / cm  CIO) 


°P*C»0+M  5=£OOT(*»)+M+102  1/cm 


C  11) 


where  *  represents  the  vibrational  energy  particle;  the  symbol 
inside  the  brackets  indicates  the  CO  vibration  mode.  The 
energy  difference  is  taken  from  translation  or  rotation. 


laser  system.  -  2  2 


Only  a  few  of  the  lower  vibrational  energy  levels  of  the 
system  are  presented  in  the  formulae  and  figures.  If  the 
relaxation  equation  of  vibrational  energy  is  used  to  describe  the 
system,  the  effects  of  all  similar  processes  can  be  included  1201; 
if  an  equation  which  considers  the  particle  number  only  at  the 
lowest  energy  level  is  used  to  describe  the  relaxation  process, 


errors  can  occur.  However,  for  the  temperature  range  (<1500K)  in 
which  we  are  interested,  this  approximation  is  acceptable. 

The  relaxation  mechanism  of  a  molecule  system  is  extremely 

complicated.  Herzfeld  C 111  has  analytically  calculated  in  detail 

for  the  CO  system  based  on  the  SSH  theory  C51  and  gives  the 
2 

Jumping  probability  for  all  possible  processes.  However,  the 
experiment  results  indicate  that  the  SSH  theory  has  a  lot  of 
constraints.  The  fact  is  that  the  jumping  probability  among 
various  vibrational  states  is  difficult  to  measure  directly,  and 
many  people  are  currently  devoting  themselves  to  determining  the 
relaxation  mechanism  indirectly  by  using  a  laser  method. 

The  .*»-**»+»%  mechanism  is  adopted  to  study  the  V-V  exchange 

inside  a  OOk  (**-** molecule  while  the  mechanism  is 

.  -  4-  1  , 

A  • 

used  irv-Ref .  C 13 .  ' 


.OOSW+M 

«=±OOr(*0+M+416 
tcM-OO*  N*  HfO) 


1/cm 


C  12) 


Under  an  equilibrium  condition  between  <XV'i)  and  modes,  it 

'  — TT1  - 

can  be  proved  that  there  exists  a  simple  relation  between  the 
rates  in  Equation  (8)  to  CIO)  and  C12)  [221, 


r  C2(0-*3),  Jf] 

x  K  [C8(l-+0),  01(0  -*1), 

C2(0-*>1),  M]  C13) 


where  0,  is  a  vibrational  characteristic  temperature.  Thus  only 


the  reaction  rates  of  Equation  (8)  to  CIO)  are  given  here. 


From  the  experiment  data  analysis,  it  shows  that  the  rate  of 


c  i4) 

"""  '  *  —  ~  .  *... ..  -  •m-v ■  •  •  &  • 

is  zero  in  the  error  range  C23],  thus  it  is  ignored. 


A  small  amount  of  water  has  a  great  effect  on  the  relaxation 

process  of  a  CO  -N  system.  Due  to  the  difficulties  of 
2  2 

experiments,  imperfect  theory  and  the  complications  of  water 

reaction,  the  relaxation  mechanism  after  adding  water  is  still  not 

clear.  In  summary,  there  are  two  viewpoints:  one  states  that  a 

*  * 

V-V  exchange  first  occurs  between  CO  and  N  ,  and  then 
.  *  '  2  2 

see  Ref.  [11);  the  other  one 

collision  companion  in  the  V-T 

*  * 

or  V-V  excitation  of  CO  ,  N  .  The  latter  theory  is 

2  2  * 
adopted  in  this  paper.  If  the  V-T  excitation  rate  for  HO  is 

2 

extremely  great,  it  can  be  proved  that  the  water  vibration  mode 
can  be  treated  as  in  equilibrium  with  the  translation,  thus  the 
effects  from  both  viewpoints  are  equivalent. 


Ill  Experiment  and  Data  Process 


There  are  many  experimental  methods  for  studying  vibrational 
relaxation.  The  most  common  methods  are  a  shock  tube 
(compressible  flow),  laser  fluorescence,  laser  photosound, 
ultrasound,  expansion  flow  in  a  shock  tube,  collision  tube, 


the  H  iQ&Jfcerfarms  V-T  excitation  C 

-  2™  M 

considers  the  Water  molecule  as  a 


6 


etc.  . 


The  relaxation  rate  is  a  microscopic  datum  which  describes  one 
of  the  characteristics  of  a  molecular  vibrational  energy  level  and 
cannot  be  measured  directly  from  experiments.  The  measurement  in 
an  experiment  is  some  sort  of  macroscopic  physical  quantity  which 
has  a  close  relation  to  the  relaxation  process  of  interior 
degree-of-f reedom.  A  certain  theoretical  analysis  is  necessary  to 
connect  the  experiment  data  Csuch  as  decay  constant  or  total 
relaxation  time)  with  the  relaxation  rate  of  some  processes. 
Obviously,  the  speed  accuracy  also  depends  on  data  process  in 
addition  to  the  experiment  itself. 

The_Schafer  formula  for  the  V-T  rate  of  a  multi-atom  molecule 
[243  l»-always  used, 


( 15) 


where  the  subscript  1  represents  the  mode  of  "control  relaxation 
rate";  C  is  the  vibrational  specific  heat.  Study  [223  shows 


& 


that  this  formula  is  good  only  when  the  V-V  rate  in  the  system  is 

much  greater  than  the  V-T  rate,  a  "control  rate";  however,  the 

actual  V-V  rate  is  not  sufficiently  large.  The  error  is  certainly 

introduced  when  the  Schafer  formula  is  adopted,  especially  under 

the  condition  of  small  specific  heat.  This  error,  however,  is 

small  in  a  CO  system. 

2 


The  V-V  rate  process  was  recently  developped.  Under  the 


f 

i 

f 

t 

I 

t 


rendition  of  a  pure  CO  and  CO  -H  0  mixture,  the  decay 

2  2  2 

constants,  k,  which  is  measured  by  a  shock  tube  experiment 


[  25  J 


performed  by  Taylor  and  Bitterman,  and  the  laser  fluorescence 
experiment  carried  by  Rosser  et.  al.  [23,  263  are  then  converted 
to  rate  coefficient  through  the  formula, 


jBT[C3(1-*0),  G2(0-»8),J(]-JL 

- - ..  26 


C  16) 


Based  on  the  common  form  of  the  relaxation  equations  [203,  a 
relaxation  mechanism  is  analyzed.  A  solution  is  also  obtained  by 
solving  a  set  of  relaxation  equations  [223, 


C  17) 


where 


#"*&&*)*  •  C19) 


(20) 


According  to  Ref.  [203,  the  factor  2  in  the  Taylor  and  Bitterman 
formula  is  combined  into  the  constant  K.  If  the  FCCC2C1  -♦  0)3  in 
Equation  (17)  is  much  greater  than  k,  it  can  be  reduced  to 
Equation  (16).  This  is  the  exact  approximation  adopted  by  Taylor, 


Bitterman  and  Rosser  which  considers  that  CO  ML  and  translation 

2 

are  in  balance.  The  rate  obtained  by  using  Equation  C 16)  will 

have  an  error  in  the  relation  of  its  numerical  value  and  the 

temperature.  It  is  more  complicated  for  the  CO  -N  mixture, 

2  2 

both  the  rate  of  C0^»r>  and  the  rate  of  ocsto-**; 

have  to  be  solve  simultaneously.  A  similar  approximation  was 
adopted  by  them,  and  thus  errors  were  created.  The  results 
processed  by  these  two  formulae  are  presented  in  this  paper. 

3 

The  unit  of  rate  K  used  in  this  paper  is  cm  /molecule- sec. 

The  relation  between  K  and  the  jumping  probability  P,  the  formula 
which  converts  various  relaxation  experiment  data  to  K,  and  those 
necessary  physical  constants  are  listed  in  Reference  Cl]. 

IV  Relaxation  Speed  Data 

Since  the  vibrational  relaxation  time  is  short,  and  the 
requirements  of  the  experiment  Csuch  as  impurity  content)  are 
strict,  the  data  measured  from  experiments  often  show  large 
errors;  some  time  these  errors  are  one  or  two  orders  of 
magnitude.  According  to  the  data  reliability  and  experiment 
accuracy  analyses,  a  certain  weight  is  added  to  various  data 
measured  by  different  authors,  and  then  a  weighted  least  square 
method  is  used  to  get  the  best  fit  in  the  form  of 


lognZ -A+BX+OJP+DX* 


(21) 


Call  calculations  are  performed  in  a  state-made  computer  TQ-16). 


According  to  Landau-Teller  and  SSH  theory,  generally,  the  X  is 
-1/3 

chosen  as  T  .  Under  a  high  temperature  condition,  should 

-1/3 

be  linear  to  T 


Table  1  Best  fitting  coefficients 


<1>  (2)  (3) 

r -  I*  *  I  r — ———I T T I  ~  "I 


** 

[7 x  ■ 

By 

0 

.  D  : 

B 

ES 

-7.22788741 

203.862916 

•  '  0 

S3 

-908704466 

. 

46.4011141 

■ 

160—1600 

’X* 

T-i^r 

-91827809 

-25.22683 

126.6784 

0 

300—1000 

(A  -N* 

»  • 

X 

-5.88873 

-189.0698 

88.50646 

1412.096 

240—4700 

WMi(r-^n-c>o5^ 

‘  X  ' 

-6.38878 

-139.0598 

88.50646 

1412.096 

240—4700 

.  (6)  Kj(F-»Tl  — HjO.  ' 

X 

-10.9764695 

-21.7181799 

0 

0 

290—3000 

X 

SE 

-12.0624762 

-0.098067132 

5.539994X10-* 

-6.277332  X  10-< 

150—2600 

v'ir 

atocL  -co. 

&Tf 

r-v» 

—9.26928286 

-4.69827096 

-470.106702 

1980:79807  , 

150—1600 

•  X 

-9.7767186 

2.92817157 

-604.993077 

1978.16692 

150—1600 

j  *  &  '  '  : 

1  «t  «>»(.*■*•«.+»*> -N, 

ESI 

r-»/» 

5.61976614 

-847.262296 

2067.64668 

-4288.96898 

160—1600 

X 

r-i/» 

-890.990028 

2387.44778 

-6029.38366 

160—1600 

Kr* 

r-v* 

-14.1877762 

88.1608647 

-159.940244 

0 

300—1000 

* 

l-v* 

-18.8868966 

38.8333884 

-144.499031 

0 

300—1000 

n 

(U>  oot^a^-oo, 

15 

UJ 

D 

-13.07283 

0.00137333 

0 

0 

300—1000 

iW5»K-i.acxio-«,l^«iMrac3tJ6)  . 

300,  400 

(U)  00,(«r-»2i^)  -N,  < 

•  •  _ 

co1(t»i-»  2v»)  —  QOj  *;**£&«*  ( 7 ), 

400 

400 

OD  00,(n-^»*)-Ho 

KIM  |9| 

Key:  ci)  Process;  C2)  Rate;  C3)  Temperature  range; 
(4)  Slow  K11C5)  Fast  K;  C6)  Considered  as 
K=l. 25x10  ,  and  not  varied  with  temperature;  C7) 

Considered  as  equal  to  the  oo,(*^* -oo*  ■  rate;  C8) 
Considered  as  0.46  time  of  -co,  rate;  C9) 

Lack  of  data. 


The  best  fit  coefficient  is  list  in  Table  1.  Figure  2  to  10 


are  plots  of  the  relation  between  rate  coefficients  and 
temperatures. 

Anderson  [18]  investigated  the  effects  of  various  relaxation 
rate  uncertainties  on  the  GDL  performance  by  using  his  simplified 
relaxation  model  and  the  base  of  the  linear  relaxation  equation. 

He  concluded  that  the  reactions  of  Equation  C3),  CIO),  (6)  and  C9) 
have  the  most  effect. 


Fig.  6  afc<n>-*».  rate 

A-[44];  +— £67]; 

A— [88];  O— [89];  #—[70];  Q— [71] 
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